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ABSTRACT 

Context. The role of magnetic fields during high-mass star formation is a matter of fierce debate, yet only a few direct probes of 
magnetic field strengths are available. 

Aims. The magnetic field is detected in a number of massive star-forming regions through polarization observations of 6.7 GHz 
methanol masers. Although these masers are the most abundant of the maser species occurring during high-mass star formation, most 
magnetic field measurements in the high-density gas currently come from OH and H2O maser observations. 

Methods. The 100-m Effelsberg telescope was used to measure the Zeeman splitting of 6.7 GHz methanol masers for the first time. 
The observations were performed on a sample of 24 bright northern maser sources. 

Results. Significant Zeeman splitting is detected in 17 of the sources with an average magnitude of 0.56 m s"' . Using the current best 
estimate of the 6.7 GHz methanol maser Zeeman splitting coefficient and a geometrical correction, this corresponds to an absolute 
magnetic field strength of 23 mG in the methanol maser region. 

Conclusions. The magnetic field is dynamically important in the dense maser regions. No clear relation is found with the available OH 
maser magnetic field measurements. The general sense of direction of the magnetic field is consistent with other Galactic magnetic 
field measurements, although a few of the masers display a change of direction between different maser features. Due to the abundance 
of methanol masers, measuring their Zeeman splitting provides the opportunity to construct a comprehensive sample of magnetic fields 
in high-mass star-forming regions. 
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1. Introduction 

. Although massive stars play an important role in the chemical 
' and energetic evolution of their host galaxies, their formation 
' mechanism remains elusive. This problem is the topic of ex- 
tensive observational and theoretical efforts. Even though few 
of the cuiTent simulations include magnetic fields, the influ- 
, ence of magnetism on the star formation processes is exten- 
■ sive as it can support a molecular cloud against collapse, af- 
' feet core fragmentation and change the feedback processes (e.g. 
Krumholz & Bonnell, 2007, and references therein). 

Most cuiTent high-mass star formation magnetic field infor- 
mation comes from H2O and OH maser polarization observa- 
tions. The observations of the H2O maser Zeeman effect using 
Very Long Baseline Interferometry (VLBI) reveal field strengths 
between 10 and 600 mG, while the linear polarization mea- 
surements reveal a complex but of t en orde red magnetic field 
morphology (e.g. IVlemmings et al.L l2006ah . Aside from H2O 
masers, tracing high-density regions (whj ~ 10** - 10" cm"-'), 
the magnetic field in the less dense surrounding regions is 
typically probed b y pola rimetric OH maser observations (e.g. 
iBartkiewicz et al.L l2005h . These observations reveal fields of 
a few mG as well as ordered structure in the magnetic field. 
However, the strongest and most abundant of the high-mass star 
formation region masers arises from the 6.7 GHz 5i - 6oA^ 
methanol transition, and for this maser hitherto only very few 
polarization observations exist. 
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Over 95% of the class II 6.7 GHz masers have been found 
to harbor warm dust emission, even though only some of the 
masers are asso ciated with a d etectable ultra-compact (UC) 
HII region (e.g. iHill et all l2005h . This seems to indicate that 
the masers probe a range of early phases of massive star for- 
mation. The 6.7 GHz masers are likely pumped by a com- 
bination of collisions and emission from nearby, warm {T > 
150 K) dust, but themselves arise from gas at much lower 
temperature (T < 50 K) with high hydrogen number densi- 
ties («Ht > 10^ ^t^T^^ ^ ^^Sh me thanol abundance (e.g. 
ISobolevet all Il997t ICragg et alLlIOOSh . Like H2O, methanol 
is a non-paramagnetic molecule, and thus both the linear and 
circular polarization fractions are small. The first polarization 
measurements were made with the Australia Telescope Compact 
Array (ATCA) on the 6.7 GHz ma ser toward a ha ndful of 
southern massive star-forming regions (lEIIingsenLl2002h and lin- 
ear polarization between few and 10% was detected. The first 
high angular resolution linear polarization maps were recently 
made using MERLIN ( Vlemmings et al., 2006b) and the Long 
Baseline AiTay (Dodson 2008). These observations indicate a 
typical linear polarization fraction of ~ 2 - 3%. In addition to the 
first linear polarization measurements, a marginal possible de- 
tection of circular polarization, caused by the Zeeman effe ct, was 
made for the masers of the ON 1 starforming region ( Gre en et all 
I200l . 

This paper presents the first significant detection of Zeeman 
splitting in the 6.7 GHz maser transition for a sample of north- 
ern massive star-forming regions. The observations, data reduc- 
tion and error analysis are discussed in §|2]and methanol maser 
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Zeeman splitting in § [3] The results are given in § |4] and are 
discussed in § |5] 

2. Observations & analysis 

2.1. Effelsberg observations & data reduction 

The 6668.512 MHz (5i -6oA^) methanol maser line of a sample 
of massive star-forming regions was observed on Nov 12th 2007 
using the 5 cm primary focus receiver of the 100-m Effelsberg 
telescope. The data were taken in position switch mode with a 
2 minute cycle time. The full width at half-maximum (FWHM) 
beam of the telescope is ~ 120" at the maser frequency. Data 
were collected using the fast Fourier transform spectrometer 
(FFTS) using two spectral windows, corresponding to the right- 
and left-circular polarizations (RCP and LCP). The spectral win- 
dows of 20 MHz were divided in 16384 spectral channels, result- 
ing in a ~ 0.055 km s 'channel spacing and were centered on the 
local standard of rest (LSR) source velocities. 

The data were reduced using the Continuum and Line 
Analysis Single-dish Software (CLASS) package. The ampli- 
tudes were calibrated using scans on 3C123, 3C286 and 3C84. 
No special efforts were undertaken to model short time-scale 
gain fluctuations during the observations, however, the small 
gain difference between the RCP and LCP was corrected. Small 
RCP and LCP gain differences do not affect the Zeeman split- 
ting determination when using the cross correlation method, as 
with this method the Zeeman splitting measurements are not re- 
liant on accurate absolute fluxes. It is found that in general the 
measured peak fluxes agree well with the literature v alues in the 
methanol maser catalogue of iPestalozzi et alj ( l2005i) . From the 
literature comparison and from a comparison of the measured 
flux of G9. 62+0. 20, which is the target of regular monitoring, 
with the flux obtained on Nov 12th with the Hartebeesthoek ra- 
dio telescope (S. Goedhart, priv. communication), the absolute 
flux errors are estimated to be ~ 10%. Still, for a few sources 
the measured flux density is up to 50% different from the pub- 
lished values, which is likely due to variability and/or structural 
changes in the maser themselves. 

2.2. The source sample 

The sampl e was taken from the 6 .7 GHz methanol maser cat- 
alogue by IPestalozzi et al.l (l2005l) and consists of some of the 
strongest (> 50 Jy) northern maser sources observable from 
Effelsberg. Priority was given to sources previously observed at 
high resolution or with existing OH/H2O maser polarization ob- 
servations. To detect the Zeeman splitting the goal was to reach 
a signal-to-noise of > 3000 and thus the total integration time 
per source was variable. As the observations were limited to 
13 hours, source selection also depended on observability and 
was occasionally adjusted based on the flux detected in a first 
observing scan. As a result, the noise level is different for each 
of the sources and ranges from 20 mJy to 100 mJy for each po- 
larization. 

2.3. Zeeman splitting determination 

The circular polarization fraction arising from the small Zeeman 
splitting of the paramagnetic methanol molecule, is extremely 

' The 100-m telescope at Effelsberg is operated by the Max-Planck- 
Institut fiir Radioastronomie (MPIfR) on behalf of the Max-Planck- 
Gesellschaft (MPG) 



low (< 0.5%). Detection of such small circular polarization frac- 
tion is very sensitive to an accurate relative calibration between 
the RCP and LCP signals. Zeeman splitting can be detected 
relatively str aightforward using the cross-correlation method, 
described by iModiaz et al.l (12005 ). This method directly deter- 
mines the Zeeman spUtting from the RCP and LCP data. It has 
been shown to be robust against the relative RCP and LCP gain 
calibration errors while still reaching a comparable sensitivity 
as performing S-curve fitting to the circular polarization spec- 
trum. Also, essential for single dish work, it is able to measure 
Zeeman splitting in the case of spectrally blended maser features 
in a straightforward manner without the need for fitting a number 
of Gaussian components to identify individual maser features. It 
was for instance used on Green Bank Telescope data to deter- 
mine limits on the Zeeman s plitting for two H2O maser me ga- 
maser galaxies dModiaz et al ., 2005; Vlemmings et al.L[2007h . 

As many of the sources will have blended maser features, 
and to reduce the effect of remaining RCP and LCP gain cal- 
ibration errors, the Zeeman splitting is thus determined using 
the cross-correlation method. Applying this method to the en- 
tire maser spectrum essentially means that the resulting Zeeman 
splitting is a flux weighted average of the true Zeeman split- 
ting. To examine changes of the magnetic field over the maser 
spectrum, one can perform the cross-correlation analysis on sub- 
sections of the spectrum. In the case of these observations, rms 
noise considerations limit the useful spectral range over which 
the analysis can be applied to > 3 km s '. Thus, to describe 
the magnetic field changes over the maser spectrum, cross- 
correlation was done in the interval [(V, - 1.5),(y,-i- 1.5)] km s"', 
with y, taken along the spectrum at steps of 3 channel widths. 
Examples of such 'running' cross-correlation are shown in Fig.[T] 
with only magnetic field values plotted that had > 5cr signifi- 
cance. Special care has to be taken with velocity intervals which 
are dominated by the wings of a bright emission peak. In this 
case, remaining gain errors between RCP and LCP could cre- 
ate a spurious Zeeman detection. Although a test with varying 
gain factors indicates that the 'running' cross-correlation is ro- 
bust against cases with a uniform gain error across the spec- 
trum, small baseline variations can still cause errors. Thus, only 
Zeeman determinations in velocity intervals that included a sep- 
arate maser feature (i.e. a local flux maximum) in the inner 80% 
of the interval were taken to be significant. It should be noted 
that, as the 'running' cross-correlation effectively convolves the 
spectra down to a 3 km s"' spectral resolution, apparent gradual 
magnetic field variations are introduced when spectral features 
with a different field strength are blended together. This grad- 
ual magnetic field change between the maser features is thus not 
necessarily physically the case in the maser region. 

Although the analysis was performed using the cross- 
correlation method, the circular polarization pattern as expected 
from Zeeman splitting was detected for several of the least com- 
plex sources. The magnitude of the circular polarization was 
~ 0.2%. Fig. |2] shows an example of two sources. The case of 
Cepheus A illustrates the problem with determining magnetic 
fields from the circular polarization produced by such a complex 
maser. Although the structure of the Cepheus A circular polar- 
ization is similar to the total power derivative, it is extremely 
difficult to properly distinguish the individual features. 

2.4. Error analysis 

Small amplitude non-Gaussian effects on the spectral baselines 
can introduce an additional uncertainty in the Zeeman splitting 
determination. To estimate this effect, Monte-Carlo modeling 
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Fig. 1. Total intensity spectra (bottom) and magnetic field strengths(top) for six sources of our sample. The magnetic field strength 
is determined from the measured Zeeman splitting using the current best value for the 6.7 GHz methanol maser splitting coefficient. 
The Zeeman splitting is derived using the 'running' cross-correlation method (see text). 



was performed using artificially generated maser spectra with 
actual baselines taken from emission free spectral regions for 
each of the sources of the sample. It was found that typically 
the rms from the cross-correlation Zeeman splitting calculations 
need to be increased by ~ 15% to accommodate this effect. This 
additional source of error has been included in the errors quoted 
below. 

2.5. Beam Squint 

The effect of a slight difference in the pointing center between 
the RCP and LCP, as a result of both feeds being off-axis, is 
called beam squint (e.g lHeilei, [T996). Beam squint can mimic 
the effect of Zeeman splitting when a velocity gradient is present 
across extended emission of a spectral line, as the RCP and LCP 
telescope beam will then be probing gas at a slightly different 
velocity. This severely complicates the Zeeman splitting obser- 
vations of extended thermal line emission and needs significant 
additional effort by rotating the feeds during the observations. 

Howe ver, typically, indivi dual maser features are not very 
extended. iMinier et al.l (|2002|) found that the majority of the 
methanol masers consist of a compact core with a diffuse 
halo structure of up to a few hundred AU, corresponding to 
up to 100 mas. Even for the most extended sources, such as 
W3(OH), where methanol maser emission extends over sev- 



eral arcseconds, the bulk of the maser flux which dominates the 
Zeeman splitting determination, comes from a region smaller 
than 100 mas. The pointing difference between the E ffelsberg 
RCP and LCP receivers are < 1" (e.g. iFiebid. Il990 l). The ef- 
fect of beam squint will be most severe when the velocity gra- 
dient exists in the direction of the beam off-set. Thus, during 
longer observation scans, any artificial circular polarization sig- 
nature will be quenched when the feeds rotate under the source. 
However, as the observation scans were occasionally as short 
as a few minutes, it was determined to what level beam squint 
could contribute to the observed Zeeman splitting. Assuming the 
worst case, a velocity gradient across the maser in the direc- 
tion of the RCP-LCP beam off-set, a broad maser spectrum of 
~ 5 km s"' and a maser extension of 100 mas, the contribution of 
the beam squint to the Zeeman splitting is found to be less than 
0.002 m s corresponding to 0.04 mG. The beam squint circu- 
lar polarization signature and the regular maser Zeeman splitting 
only becomes similar if the maser would extend over more than 
5". However, when an extended maser region consists of indi- 
vidual compact maser features that dominate the spectrum, the 
beam squint effect is negligible. It can thus be concluded that 
the velocity splitting between RCP and LCP measured for the 
6.7 GHz methanol maser sources is not due to beam squint. 
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Fig. 2. Total intensity and circular polarization spectrum for G37.40+1.52 (left) and G109.86+2.10 (Cepheus A; right). The thick 
solid line in the bottom panel is best fit fractional total power derivative to the circular polarization spectrum 



3. Methanol maser circular polarization 

3.1. Zeeman splitting 

The methanol molecule is non-paramagnetic and as a result the 
Zeeman splitting under the influence of a magnetic field is ex- 
tremely small. The split energy, AEz, of an energy level un- 
der the influence of a magnetic field, B, is described by AE, - 
IQ'^gfiffMjB. Here My denotes the magnetic quantum number 
for the rotational transition with the rotational quantum number 
J, B is the magnetic field strength in Gauss, fiN is the nuclear 
magneton and g is the Lande g-factor, which thus determines 
the magnitude of th e Ze eman eff'ect. This factor was investigated 
many years ago by IJenI (Il95ll) . who found empirically that it is 
probably an average of the true g-factor of several interacting 
states and can be described by the equation: 

^ = 0.078 + 1.88/[7(7h-1)]. (1) 

This is based on laboratory measurements using 25 GHz 
methanol masers lines and as a result one has to be careful apply- 
ing these results to the 6.7 or 12.2 GHz masers. However, until 
further laboratory studies, this remains the best estimate and im- 
plies, for the 6.7 GHz methanol transition, a Zeeman splitting 
coefficient of 0.049 km s"' G"'. It is difficult to properly judge 
the uncertainty in this value. Conservatively, one expects that 
magnetic field strengths determined using this value are within 
~ 25% of the true value. 

3.2. Non-Zeeman effects 

There are several, non-instrumental, effects other than Zeeman 
splitting that might cause circular polarization of the maser line. 
One possibility is a rotation of the axis of symmetry for the 
molecular quantum states. This can occur when, as the maser sat- 



urates, the rate for maser stimulated emission R becomes larger 
than the Zeeman frequency shift gQ.. While gQ. » R, the mag- 
netic field direction is the quantization axis. Then, when R be- 
comes smaller than gQ., the molecules interact more strongly 
with the radiation field than with the magnetic field and the 
quantization axis changes towards the maser propagation di- 
rection. From the Zeeman splitting coefficient derived above, 
gQ. X 22Z?[mG] s ' for the 6.7 GHz methanol maser The rate 
for stimulated emission can be estimated using: 

R AkTbAQ/4nhv. (2) 

Here A is the Einstein coeffic ient for the maser t ransition, which 
is equal to 0.1532 x 10 s"' (Cragg et al.L ll993h . and k and h are 
the Boltzmann and Planck constants respectively. The maser fre- 
quency is denoted by y, and T\, and AQ are the maser brightness 
temperature and beaming sohd an gle. Observations indicate that 
Tb < 10'^ K (Mi nieret aLll2002h . but AQ is somewhat harder 
to estimate and decreases rapidly with increasing maser satura- 
tion level. For a typical maser beaming angle of AQ » 10"^, the 
maser stimulated emission R < 4 s '. Thus, even in the most 
saturated case gQ » R when the magnetic field strength is 
> 1 mG. A change of quantization axis is thus unable to explain 
the observed circu lar polarization. 

Alternatively, Wiebe & WatsonI (Il998l) have shown that the 
propagation of strong linear polarization can cause circular po- 
larization when the direction of the magnetic field changes sig- 
nificantly along the maser propagation direction. For a smooth 
change of magnetic field direction of ~ 1 rad along the maser, the 
fractional circular polarization caused by this effect is approxi- 
mately mf/4, where mi is the fractional linear polarization. For 
the typical 6.7 GHz methanol maser linear polarization fraction 
of 2 - 3%, this implies that the generated circular polarization is 
only ~ 0.02%, which is less than 10% of the observed circular 
polarization. And, the fairly constant linear polarization vectors 
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observed at high-resolution (e.g. lVlemmings et an,l2006bl) . indi- 
cate that the magnetic field rotation along the maser path length 
is likely much smaller than 1 rad, making the propagation of 
linear polarization an unlikely source of the methanol maser cir- 
cular polarization. 

Finally, a velocity gradient in intrinsic circular polarization 
(i.e. not caused by the magnetic field) could mimic the effect of 
Zeeman splitting. However, as in a maser the stimulated emis- 
sion quickly dominates the intrinsic emission, any intrinsic cir- 
cular polarization qu ickly becomes negligible. This is described 
in more detail in e.g. lElitzuii (Il998h . 

4. Results 

The results of the Zeeman splitting analysis using the RCP-LCP 
cross-correlation method are presented in Table[T] The table lists 
the source name, position, central Vlsr velocity, peak and inte- 
grated flux and the measured Zeeman splitting AV^. As described 
in § 12.31 this value corresponds to the flux averaged Zeeman 
splitting of the entire maser spectrum. In the final column it fists 
the strength of the magnetic field component along the line-of- 
sight determined using the current best value of the 6.7 GHz 
methanol maser Zeeman splitting coefficient as discussed in §|3] 
Throughout the paper, the quoted magnetic field strength will 
correspond to this value. Should in the future a better value of the 
Zeeman coeflicient be determined then all magnetic field values 
need to be adjusted correspondingly. 

Fig. [T] shows six methanol maser spectra and correspond- 
ing magnetic fields derived using the running cross-correlation 
method described in § 12.31 The figure shows that for a number 
of the sources, the magnetic field changes across the maser. For 
those sources. Table [T] also lists separately the magnetic field 
derived for sub-sections of the spectrum. Total intensity spec- 
tra and magnetic field strengths for all the sources in the sample 
with a significant Zeeman detection are presented in the Online 
material as Figs.|6]and|7] Online Fig.[8]shows the total intensity 
spectra for the sources without Zeeman detection. 

4.1. Individual Sources 

A number of the maser sources from the sample have been ob- 
served at high angular resolution, providing information on the 
morphology of the methanol maser region. Additionally, several 
sources have been the target of magnetic field measurements 
using mainly OH and/or H2O maser polarization observations. 
Sources for which high-resolution 6.7 GHz observations and/or 
additional magnetic field measurements are available are dis- 
cussed here in detail. 

4.1.1. 09.62+0.20 

The 6.7 and 12.2 GHz methanol masers associated with the high- 
mass star formation complex G09.62-h0.20 undergo periodic 
flares with a period of 244 days, making these the first reported 
incidence of periodic variations in a high-mass star-forming re- 
gion (Goedhart et al., 2003, 2004). The origin of the periodic be- 
havior, however, is still unclear. The 6.7 GHz methanol m asers 
have been mapped with the ATCA bv lPhiUips et all (11998 ). and 
the strongest features are shown to be associated wi th the hyper- 
compact HII region labeled E by Garav et aTl ( Il993h . This region 
also hosts OH masers for which Zeeman sph tting observation s 
indicates a magnetic field strength of ~ 5 mG dFish et alll2005[) . 
As seen in Fig[T] the magnetic field determined from the 6.7 GHz 



methanol maser Zeeman splitting reverses direction between the 
second-strongest maser feature at Visr = -0.1 km s"'(B|| = 
8.6 mG) and the strongest feature at 1.2kms"'(B|| - -2.96 mG), 
even though the ATCA image shows that the features are located 
within ~ 100 mas of each other. The periodicity of the flaring 
of both maser features is the same, with the weaker of the two 
features flaring with a ~ 1 month delay. There is an intriguing 
possibility that the maser flaring and time delay is the result of 
a periodic change in the ma gnetic field or i entatio n due to a form 
of magnetic beaming (e.g. Grav & Field |1994|) . However, this 
hypothesis needs to be tested by magnetic field monitoring ob- 
servations and improved maser modeling. 

4.1.2. 12.89+0.49, IRAS 18089-1732 

This young massive star-forming region contains methanol, 
H2O and OH maser emission as well as a weak hypercom- 
pact HII region. The 6.7 GHz methanol masers have been 
mapped with the ATCA and c onsist of two regions separated 
by ~ 1.5"(\^lihetal], Il997h . The strong maser feature at 
Vlsr = 33.6 km s"', which shows the smallest Zeeman split- 
ting (Z?|| - 4 mG) is l ocated closest to t he sub- millimeter contin- 
uum peak imaged bv'B euther et al.l (e.g. l200 5h at high-resolution 
with the sub millimeter array (SMA). The higher velocity fea- 
tures, with B|| - 12 mG, are located in the direction of the 
North-South molecular outflow found in the same SMA obser- 
vations. The OH masers on the other hand, are located in the 
direction of extended molecular emission in the East- West direc- 
tion, which shows indications of a rotat ion signature perpendicu - 
lar to the outflow dBeuther et all 120051: iBeuther & Walslil2008l) . 
The magnetic fi eld measured from OH ma ser Zeeman splitting 
isB = 3.93 mG dRuiz-Velasco et al.Ll2006l) . 

4.1.3. 23.01-0.41 

The Zeeman splitting for this source, and consequently the mag- 
netic field, shows a very strong negative magnetic field for the 
brightest maser peak at Vlsr - 74.9 km s"'(B|| - -30.3 mG), 
while the other maser features have a positive field direction with 
a strength of By ^ 14 mG. In Fig. [1] this results in a dip to- 
wards an increasingly negative magnetic field across the maser, 
because, as described in § 12.31 the running average calculates the 
Zeeman splitting over a 3 km s~ 'interval. ISzvmczak & GerardI 
(2004') find a 1667 MHz OH maser feature with a magnetic field 
strength of B = -0.2 mG at Vlsr = 74.3 1 km s"' . Unfortunately, 
no high resolution maps are available to check the morphology 
of the maser source. 

4.1.4. 31.28+0.06, IRAS 18456-0129 

The largest Zeeman splitting was measured for the masers of the 
UC HII region G3 1 .28+006, which is part of the giant HII region 
W43, and correspond to By - 42 mG. The 6.7 GHz methanol 
masers were observed with the European VLBI Network (EVN) 
and reveal a complex distribution over ~ 400 x 400 mas 
dMinier et all l200d) . No other maser polarization observations 
are found in the literature. 

4.1.5. 33.64-0.21 

As seen in Fig. [1] the magnetic field of this source is stable 
over the entire velocity range with a flux averaged strength 
of B|| - -18 mG. The only other magnetic field measure- 
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Table 1. Zeeman splitting results 
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IRAS 18517+0437 


18 54 13.80 


+04 41 32.00 


41.0 


320 


193 


0.75 ± 0.02 


15.4 ±0.4 


43.80-0.13 


W49N 


19 11 55.10 


+09 36 00.00 


40.0 


35 


81 


< 1.9 


< 39 


49.49-0.39 


W51-el/e2 


19 23 44.50 


+ 14 30 31.00 


59.0 


1029 


885 


0.72 ± 0.04 


14.7 ±0.8 


69.52-0.97 


ONI 


20 10 09.07 


+3131 34.40 


11.6 


96 


45 


<0.2 


< 4.1 


78.10+3.64 


IRAS 20126+4104 


20 14 26.04 


+41 13 33.39 


-6.1 


60 


77 


< 0.8 


< 17 


81.87+0.78° 


W75N 


20 38 36.80 


+42 37 59.00 


5.0 


273 


317 


0.40 ± 0.02 


8.2 ±0.3 










5.0* 






0.28 ± 0.02 


5.7 ±0.4 










7.3* 






0.46 ± 0.02 


9.4 ± 0.4 


109.86+2.10 


Cepheus A 


22 56 18.09 


+62 01 49.45 


-4.2 


364 


484 


0.39 ±0.01 


8.1 ±0.2 


111.53+0.76° 


NGC 7538 


23 13 45.36 


+61 28 10.55 


-56.1 


233 


514 


0.79 ± 0.03 


16.2 ± 0.6 










-58* 






0.78 ± 0.04 


16.0 ±0.8 










-56* 






0.59 ± 0.02 


12.1 ±0.4 


133.94+1.04 


W3(OH) 


02 27 04.72 


+61 52 24.73 


-44.0 


3705 


8198 


0.141 ±0.003 


2.87 ± 0.05 


173.49+2.42 


823 1 


05 39 13.06 


+35 45 51.29 


-13.0 


52 


53 


0.95 ±0.11 


19 ±2 


174.19-0.09 


AFGL5142 


05 30 42.00 


+33 47 14.00 


2.1 


55 


34 


< 0.8 


< 18 


188.95+0.89 


IRAS 06058+2138 


06 08 53.35 


+21 38 28.67 


10.9 


633 


485 


-0.49 ± 0.02 


-9.9 ± 0.4 


192.60-0.05 


S255 


06 12 54.02 


+ 17 59 23.00 


5.0 


94 


78 


0.47 ± 0.05 


9.6 ±0.9 



° Significant magnetic field changes across the maser spectrum. 

* Velocity of the maser feature for which Zeeman splitting was determined separately. 



ment for 33.64-0.21 comes from a single 1720 mG OH maser 
feature at Vlsr - 60.23 km s which has B - -1 mG 
jSzvmczak & Gerard l2004l) . 



4.1.6. 35.20-0.74, IRAS 18556+0136 



The methanol masers of the bipolar outflow source G35.20- 
074N hav e only recentl y been mapped using the Japanese VLBI 
Network dSugivama et al.. 2007), showing that the brightest fea- 
tures make up two compact regions separated by more than 2"in 
a dkection per pendicular to the CO outflow, along a putative 
molecular disc dPent et al.L Il985h . The Zeeman splitting obser- 
vations indicate a stable flux averaged magnetic field strength 
of B|| - 16.5 mG, which is likely dominated by the strongest 
maser feature at Vlsr = 28.7 km s^'in the southern maser clus- 
ter The OH mase r polarization of this source ha s been studied 
with MERLIN by Hutawarakorn & CohenI (Il999l) . who find that 
the magnetic field reverses on opposite sides of the disc. The OH 
maser distribution is thought to lie along the disc, with the south- 
ern masers tracing a mean field B x 4- mG and a northern maser 
feature having B = -2.5 mG. 



4.1.7. 35.20-1.74, W48 

The UC HII region G 35.20-1.74, in the Galactic HII region 
W48, is believed to be a site of massive star formation and 
has both OH and methanol masers. EVN observations of the 
6.7 GHz methanol maser region shows a ring-like structure 
of ~ 200 X 400 mas (Minieret al., 2000). Although the spec- 
trum is complex, the Zeeman splitting observations do not show 
any large variations and yield a flux averaged field strength of 
B|l - 6.4 mG. The magnetic field at the periphery of the UC 
HII region has been determined from carbon recombination line 
observations, which give a field of B = 2.9 mG at a hyd rogen 
number density of =8.5 10'' cm"^ (iRoshi et al.Ll2005l) . 

4.1.8. 49.49-0.39, W51 -el /e2 

W5 1 is one of the most luminous massive star formation com- 
plexes of our Galaxy and can be divided in three regions, W5 1 A, 
W51B and W51C OCarpenter & Sandersl[l998l) . The maser site 
G49.49-0.3 9 is associated with W51A and particularly with 
W51-el/e2 jCaswell et al.Lll995h . OH maser polarization obser- 
vations reveal a total of 46 Zeeman pairs near el and e2, with the 
region having predominantly a positive ma gnetic field d i rectio n 
which reverses in the northern part of el (iFish & Reidl l2006l) . 
W51e2 has two Zeeman pairs that imply a magnetic field of the 
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order of ~ 20 mG dArgon et al.L l2002h . with the rest of the OH 
masers indicating B ~ 5 mG. For the 6.7 GHz methanol masers 
the Zeeman splitting indicates a flux averaged By - 14.7 mG. 

4.1.9. 69.52-0.97, ONI 

The 6.7 GHz methanol maser emis sion and polarizat ion was 
recently mapped with MERLIN by iGreen et all ( l2007h . Those 
observations indicate a tentative first detection of 6.7 GHz 
methanol maser Zeeman splitting of 9 + 3 m s correspond- 
ing to a field strength of -18 mG. This marginal detection could 
not be confirmed in the observations presented here. However, 
the Zeeman splitting measured by Green et al. was found on a 
maser feature that, in lower resolution observations, would be 
blended both positionally and spectroscopically with the bright- 
est maser. As a results, the flux averaged Zeeman splitting mea- 
surement is biased toward a possibly negligible magnetic field 
of the brightest feature. 



4.1.10. 81.87+0.78, W75N 

The 6.7 GHz methanol masers of the very active region of mas- 
sive star formation W75N have been mapped with the EVN 
jMinier et al.L l2000h . Their map reveals that the masers make 
up two distinct regions, an elongated region of ~ 200 mas with 
Vlsr < 6 km s^'and a compact feature ~ 500 mas South-East 
at Vlsr = 6.8 km s"^ The observations shown in Fig. [1] high- 
light the distinct nature of these maser features as the masers in 
the linearly extended structure have By =5.7 mG while those in 
the compact region have By - 9.5 mG. W75N has been the tar- 
get of numerous OH maser polarization observations indicating 




4 6 8 

log Hj, [cm-3] 



12 



another possible disc related field re versal and typical magnetic 
field streng t hs of I B| ^5-7 mG (e.g. lHutawarakorn et al.Ll2002l: 

ISlvsh et al.L 120021) . Magnetic field measurements during a OH maser measurements (solid dots; 'Vlemmings et al.', '2006 j|), av 



Fig. 3. Magnetic field strength B in the massive star forming 
region Cepheus A measured from Zeeman measurements as a 
function of riu^, the number density of neutral hydrogen. The 
solid triangle indicates the methanol maser measurement from 
the Eff'elsberg observations, the number density is taken to be 
in the range typical of 6.7 GHz methanol masers (Crag g et all 
2005). but cannot be determined precisely without observa- 
tions of additional transition lines. Also indicated are the H2O 



maser flare in W75N re veal the strongest OH m aser magnetic 
field B = 40 mG to date dSlvsh & Migenesll2006h . 



4.1 .1 1 . 1 09.86+2.1 0, Cepheus A 

The masers of the Cepheus A mass ive star formation region 
has been studied in great detail (e.g. IVlemmings et al.L l2006ai 
and references therein). The 6.7 GHz methanol maser distri- 
bution has been mapped with the Japanese VLBI Network 
jSugivama et al.L l2007i) and the EVN (Torstensson et al., in 
prep.). These observations show that the masers are found in 
an elongated structure of ~ 1500 mas across Cepheus A HW2. 
The measured methanol maser magnetic field is stable at By - 
8.1 mG. 

Cepheus A is the source in the sample for which the most 
magnetic field measurements at different hydrogen density are 
available. Fig.[3]presents on overview of all these measurements. 
It seems clear that the magnetic field strength B in Cepheus A 
has a power-law dependence on the hydrogen number density 
nH2, with the best-fit giving B oc n-O-'iT^om^ This is consistent 
with the empirical relation B oc n^"^ from low-density molec- 
ular cloud Zeeman observations by Crutcher ( 1999). Of course, 
one has to be cautious when relating the maser observations with 
the other Zeeman splitting results and the dust polarization ob- 
servations at very different scales. Additionally, the methanol 
and H2O maser magnetic field strengths depend on an assump- 
tion with regard to the angle between the maser propagation di- 
rection and the magnetic field and their number density has to be 
determined using maser excitation models. 



erage OH maser measurement (solid square; Bartkiewicz et aLl 
l2005b . NH3 measurement (open square; Garay et al., 199a) 
and recent dust polarization measurement (open triangle; 
ICurran & Chrysostomoui l2007h . The error bars on the density 
of these data are taken from the respective publications. The 
solid line is the theoretical relation B oc fixed to the OH 
maser measurement. The dashed line is the error weighted best 
fit line to the presented measurements (the individual H2O maser 
measurements were averaged). Additionally, the boxes indicate 
the range of literature values for Hl/non-masing OH, OH maser, 
methanol maser (this paper) and H2O maser Zeeman splitting 
observations for other massive star formation regions. 



4.1.12. 111.53+0.76, NGC 7538 

The methanol masers of NGC 7538 have been propo sed to 
trace a disk around a high-mass protostar (Pestalozzi et al.L 
200 4|), however, this inter pretation has been questioned by 
De Buizer & Minied (l2005h who find that the maser might be 
related to an outflow. The observations presented here indi- 
cate a flux averaged magnetic field By = 16.2 mG. The 'run- 
ning' cross-correlation however, reveals a more complex pic- 
ture. The spectrum is made up of several strong maser features 
with apparently different magnetic field strengths. The strongest 
magnetic field ~ 16 mG is measured on the maser feature at 
Vlsr = -58 km s"', while ~ 12 mG is found for the maser at 
Vlsr = -56 km s"' . In both cases the 'running' cross-correlation 
derived magnetic field is somewhat supressed as the masers be- 
tween Vlsr - -58 and -56 km s"' have a field of only a 
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few mG. Meanwhile, no significant magnetic field is detected 
in the masers near Vlsr = -61 km s"', implying a field strength 
B|l < 5 mG. The flux averaged magnetic field is marginally larger 
than the field determined using the 'running' cross-correlation 
for individual maser features. This is due to the fact that, in de- 
termining the flux averaged field strength, strong fields, which 
have nevertheless < 5cr significance in the smaller 3 km s ' in- 
tervals used for the 'running' cross-correlation, still contribute. 
OH maser measurements in dicate the region un dergoes a field 
reversal and has |B| = 1 mG (iFish & Reidll2006l) . 



4.1.13. 133.94+1.04, W3(0H) 

Linear polariz ation of the 6.7 GHz metha nol masers of W3(OH) 
is described in lVlemmings et all ('2006b'), who also give a com- 
pilation of previous OH maser magnetic field strength measure- 
ments. The polarization observations show that the magnetic 
field traces the extended methanol filament but has a more com- 
plex structure in the do minating compact region described in 
iHarvev-Smith & CohenI ( 1200 6). The Zeeman splitting presented 
here corresponds to a flux averaged magnetic field strength of 
B|l = 2.87 mG. However, due to the rich spectrum and large ex- 
tent of W3(OH) a direct comparison with other Zeeman splitting 
measurements is impossible and would need high resolution ob- 
servations. 



5. Discussion 

5.1. Magnetic field strength 

Significant Zeeman splitting was detected for 17 out of 24 
sources and indicates an absolute magnetic field strength com- 
ponent along the maser propagation direction between 2.8 
and 42 mG. There are several effects that bias the observed field 
strength towards higher or lower values. As is the case for H2O 
masers, low spatial resolution observations, velocity gradients 
across the maser and increased maser saturat ion tend to cause an 
underestimate of the magnetic field strength (Sarma et al I I2OOII 
[Vlemmings, 2006). As discussed above, beam squint can create 
a false Zeeman splitting signature, however, for milliarcsecond- 
scale masers this is unlikely to contribute more than ~ 0.04 mG. 
Finally, the Zeeman coefficient which is used to determine the 
magnetic field is uncertain, as is described in § |3] which leads 
to a systematic bias in the magnetic field. Thus the uncertainty 
in the measured absolute magnetic field strength will be sub- 
stantial, with most of the described effects biasing it towards a 
value that is lower than the actual field strength. Weighing the 
field strength by measurement significance, the average mag- 
netic field (B||,meth) = 12 mG. This needs to be corrected for 
a random angle between magnetic field and the line-of-sight, 
which implies for the absolute field strength \B\ = 2{B\\} (e.g. 
ICru tcher, 1999). This gives |Bmethl = 23 + 6 mG. The error on 
the absolute field strength is dominated by the estimated uncer- 
tainty in the Zeeman coefficient. This field strength is larger than 
that found in the OH maser regions, with t he average OH m aser 
determined field strength being ~ 4 mG ( Fish & Reidll2006l) . 

The dynamical importance of the magnetic field can be 
quantified by defining a critical magnetic field strength B^n - 
(Snpv^y^^ for which the dynamic and magnetic pressure are 
equal. Here p and v are the densi t y and velocity of the maser 
medium respectively. Crag g et"aLl (l2005l) find that the current 
observational limits suggest that the majority of the 6.7 GHz 
maser sources occur near the high-density limit of the maser 
range (ne, = 10^-10^ cm"-'), with a typical value of hh, - 
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-30 



• OH 1.6 GHz 

A OH 6 GHz 



10 

<D„„> [mG] 



20 



Fig. 4. Measured flux averaged 6.7 GHz methanol maser mag- 
netic field strength vs. average 1.6 and 6 GHz OH maser mag- 
netic field strengths for those sources where both are avail- 
able. Most literature values were taken from the compilation of 
molecular cloud magnetic field measurements bv lHan & ZhangI 
(120071) . 



10 cm . Taking this density and a typical gas v elocity of 



' from proper motion measurements (e.g. IXu et al.L 
s; 12 mG. The measured magnetic field vales are 



< 5 km s 
2006), Bcrit 

thus comparable to Bciit and hence dynamically important. 



5.2. Comparison witli OH masers 

Theoretical modeling indicates that, given sufficient abundance, 
both 1.6 and 6 GHz OH masers as well as the methanol masers 
can be pumped simultaneously (Cragg etal., 2002). However, 
observations have shown that, although they are closely associ- 
ated, co^propagation between OH and methanol masers is rare 
(e.g. lEtoka et al.. 2005). Still, as the masers are closely related, 
one might expect a relation between the OH magnetic field 
strengths and those derived from the methanol masers. As seen 
in Fig.m no clear relation is found for those sources which have 
magnetic field strength measurements from both maser species. 
This is likely due to the fact that the flux averaged methanol 
maser magnetic field strengths cannot be easily related to the OH 
maser measurements often performed at much higher resolution. 

Fig. ID does, however, show that the magnetic field direc- 
tions determined from methanol and OH masers are fully consis- 
tent. As discussed for some of the individual sources, this even 
holds for regions where the OH masers show a magnetic field 
reversal. This implies that the magnetic field orientation derived 
from the methanol maser regions, as that from OH masers, could 
be indicative of the Galactic magnetic field. This is illustrated 
in Fig. |5j although a much larger number of measurements as 
well as accurate distances to the star-forming regions would be 
needed. Since several efforts are underway to determine accu- 
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X [kpc] 



Fig. 5. Magnetic field direction and strength along the line-of- 
sight derived from the methanol maser Zeeman splitting obser- 
vations presented in this paper projected onto the Galactic plane. 
The solid dots are th e observed star-forming region with kine- 
matic distances from iPestalozzi et al.l (|2005|) These are known 
to be highly uncertain. The length of the vectors is scaled with 
B\\. The ap proximate locatio n of the spkal arms is indicated as 
taken from lTavlor & CordesI (Il993l) . 



rate distances to methanol maser regions (e.g |Xuetal.Ll2006h . 
methanol maser Zeeman splitting measurements will provide an 
important opportun ity to study the Galactic magnetic field (e.g. 
iHan & Zhangll2007l) . 

5.3. The magnetic field vs. density relation 

Fig. |3] shows the B-field measurements of both maser and non- 
maser Zeeman observations as a function of hydrogen number 
density. The figure seems to indicate tha t the .B-field follows the 
B oc n°-^ de nsity scaling law discussed in lMouschovias & Ciolekl 
over an enormous range of densities. This makes 
the maser measureme nts co nsistent with low-density molecular 
cloud measurements (ICrutc her. 1999), implying that the mag- 
netic field remains partly coupled to the gas up to the highest 
number densities. However, some care has to be taken in inter- 
preting this relation, as for instance the shock excited masers are 
short-lived (H2O masers have a typical lifetime t,„ ~ 10^ s) com- 
pared to the typical ambipolar diffusion time-scale at the highest 
densities (Td ~ 10^ s). This suggests that in the non-masing gas 
of similar densities, magnetic field strengths are likely lower due 
to ambipolar diffusion. Still, the maser magnetic field measure- 
ments strongly imply a dynamical importance of magnetic fields 
during the high-mass star formation process, with the methanol 
maser observations presented here filling the gap between the 
OH and H2O maser observations. 



6. Concluding remarks 

This paper presents the first significant Zeeman splitting mea- 
surements obtained on the 6.7 GHz methanol maser As this 
highly abundant maser uniquely pinpoints massive star forma- 
tion, a detection of the Zeeman splitting gives the opportunity 
to measure the magnetic field in a large number of high-mass 
star-forming regions at densities of ~ 10^ cm"^. The av- 
erage line-of-sight magnetic field in the methanol maser region 
(fi|i,meth) -12 mG, although this depends on the exact Zeeman 
coeflicient used to calculate the field strength as discussed in 
§ [3] A statistical correction for a randomly oriented magnetic 
field gives |Bmethl - 23 mG. This indicates that the magnetic 
field is dynamically important. The 100-m Effelsberg telescope 
observations presented here, detected significant magnetic fields 
in 70% of the sources with peak fluxes down to ~ 50 Jy. Thus, 
with additional observations it will be possible to construct a cat- 
alogue of magnetic field measurements for over 100 high-mass 
star-forming regions giving unique insight in Galactic magnetic 
fields. 
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Fig. 6. Total intensity spectra (bottom) and magnetic field strength (top) for all the sources of our sample with significant Zeeman 
splitting detection. The magnetic field strength is determined from the measured Zeeman splitting using the current best value for 
the 6.7 GHz methanol maser splitting coefficient. The Zeeman splitting is derived using the 'running' cross-correlation method (see 

§E3. 
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Fig. 8. Total intensity spectra of the sources in our sample for which no significant Zeeman splitting was detected. 



